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SUMMARY 


For  many  years,  it  has  been  recognized  that  the  amplitude  of  the  normal 
deformation  at  an  appropriate  point  on  a  column  or  on  a  plate  could  he 
connected  wj the  initial  deformation  at  that  point  by  a  simple  rela¬ 
tionship  between  actual  load  and  classic  load  for  the  structure .  The 
appropriate  expression  is  , 


6 


6 
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This  formula  is  of  considerable  importance  in  the  interpretation  of  test 
data.  However,  it  applies  only  for  small  displacements  and  requires 
that,  in  many  cases,  the  initial  imperfection  be  small  in  relation  to  the 
relevant  structure  parameter.  This  is  a  restriction,  since  it  implies 
the  need  for  high-quality  test  vehicles.  A  large  displacement  formula 
which  does  not  require  this  constraint  can  be  developed  as  follows: 


P  =  P  (1  +  y62) 
cr  '  '  ' 


This  applies,  of  course,  only  when  the  deflections  are  large  and  when 
the  load  to  produce  them  is  in  excess  of  ! -v  classic  critical  load.  A 
wide  range  of  stability  problems  is  dea7  &  v ith  in  this  report,  and  criti¬ 
cal  load  values  are  determined  from  the  appropriate  test  result  by  both 
formulae.  The  results  are  in  excellent  agreement. 
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INTRODUCTION 


The  Interpretation  of  test  data  obtained  from  instability  studies  is 
frequently  regarded  as  a  difficult  problem.  Actually,  a  powerful  method 
exists  for  its  solution.  The  method  was  first  given  for  the  strut, 2 
but  subsequent  research  has  indicated  its  generality ^  According  to 
this  method,  the  deflection  at  an  appropriate  point  can  be  related  to 
the  classic  and  actual  loads  by  the  equation 


This  is  the  equation  of  a  rectangular  hyperbola;  when  plotted  with  6/P 
and  6  as  variables,  it  gives  a  straight  line  (a  Southwell  Plot)2  whose 
slope  is  the  classic  load.  The  technique  is  restricted  to  the  inter¬ 
pretation  of  observations  made  over  a  limited  range  of  the  load-dis¬ 
placement  curve;  in  particular,  it  is  restricted  to  displacement  corres¬ 
ponding  to  loads  below  the  critical  value.  Thus,  it  is  confined  to 
behavior  within  the  limitations  of  small  displacement  theories. 

In  an  attempt  to  describe  instability  behaviour  more  adequately,  analysts 
have  developed  large  deflection  theories.  This  report  is  concerned  with 
interpretations  made  on  the  basis  of  such  treatments.  It  is  found,  as 
will  be  demonstrated  herein,  that  there  is  a  general  relationship,  applica¬ 
ble  in  the  postbuckling  range,  associating  deflection  with  critical  and 
actual  load  parameters.  Because  of  the  well-founded  analysis  of  the 
buckling  and  postbuckling  of  a  column,  the  demonstration  starts  from  this 
point. 
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COLUMN  STRUCTURES 


The  load-displacement  relationship  for  a  centrally  compressed  column  is 
shown  in  Figure  1.  In  this  diagram,  the  actual  and  ideal  motions  are 
portrayed.  The  theoretical  large  displacement  curve  is  approximated  by 
the  actual  load-displacement  relationship  when  the  deflections  become 
large.  Now,  analytically,  the  relation  between  the  load  and  the  normal 
displacement  of  the  midpoint  of  the  columns  is  expressed  by  the  following 
equation? 

.  -  M:  rooi2  b 

L2  (2) 

K 

=  f^kT  (3) 
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where  K  =  sin  ^0  and  F(K)  is  the  appropriate  elliptic  integral  of  the 
first  type  -  0  being  the  end  slope .  However,  provided  0  is  less  than  25°, 
the  complex  relationship  given  by  equations  (2)  and  (3)  can  be  well 
approximated  by  the  following  parabola :  ^ 


where 


P  =  Pcr  (1  ♦  Y62) 


(4) 

(5) 


This  parabola  can  be  plotted  as  a  straight  line  if  the  variables  are 
taken  to  be  P  and  6^.  This  line  then  cuts  the  load  axis  at  a  point 
which  corresponds  to  the  critical  value. 
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Mr.  E.  Way  carried  out  a  test  on  a  slender  column  to  verify  this  result  • 
His  data  are  given  in  Figures  2  through  4.  The  critical  load  deter¬ 
mined  from  the  slope  of  the  t/P  -  6  line  of  Figure  3  is  in  excellent 
agreement  with  the  critical  value  determined  by  the  intercept  of  the 
P  -  6*"  line  with  the  load  axis  seen  in  Figure  4.  Both  values  are  in 
perfect  accord  with  the  classic  value  computed  from  the  fornula 


P 


cr 


(6) 


A  similar  analysis  made  upon  a  more  robust  column  of  T  cross  section* 
gives  an  identical  result.  This  is  shown  in  Figure  5»  It  is  interesting 
to  note  that  this  column  experienced  some  local  yielding,  which  is 
clearly  evidenced  in  the  Southwell  Plot,  ?  Figure  6.  Nevertheless,  the 
agreement  is  still  excellent  -  7900  pounds  from  the  Southwell  Plot  as 
opposed  to  8050  pounds  from  the  6*  plot,  Figure  7. 

The  case  of  a  single  member  in  axial  compression  can  readily  be  extended 
#Test  data  from  Hill,  Reference  6. 
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Displacement 
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Figure  1.  Column  Data  with  Elastic  a  Superposed. 
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Asl&l  Load  C* 
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Figure  1*. 


Plot  of  P  vereuB  6  for  Column  Data  of  Figure  2. 
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4,  DfflKtlM  (il.) 


Figure  5.  Test  Data  from  Hill,  an  Aluminum  Column 
Buckled  Veil  Into  Fostbuckllng  Range, 
Reference  6. 
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6,  Deflection  (ia.) 
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Figure  6.  Southwell  Plot  for  Coluna  Data  of  Figure  5. 
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to  that  of  the  structure  composed  of  such  members.  For  the  Southwell-  q 
type  Plots,  this  was  established  by  the  general  analysis  of  Westergaard 
and  by  the  more  specific  analysis  and  tests  of  Gregory?*^  Recently, 
Britvec  and  Chilver^1  and  Roorda^2  £ave  published  papers  which  indicate 
the  full  applicability  of  the  P  -  6  process  for  this  class  of  structure. 
Their  demonstrations  are  based  upon  generalized  analysis  and  experiment. 
They  give  many  examples  and  deal  with  a  wide  range  of  structural  configu¬ 
rations.  The  example  we^have  chosen  to  illustrate  this  point  is  taken 
from  the  work  of  Roorda^and  is  for  the  case  of  a  two-bar  frame.  The 
details  of  the  loading  aad  the  load-displacement  relationships  for  clock¬ 
wise  and  counterclockwise  instabilities  are  given  in  Figure  8.  When  the 
appropriate  6/P  versus  6  and  P  -  6^  curves,  Figures  9  and  10,  are  con¬ 
structed,  it  is  seen  that  excellent  agreement  is  reached. 
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PLATE  STRUCTURES 


Although  the  elastica  curve  has  been  well-known  foremost  of  this  century, 
it  must  he  recorded  that  the  first  use  of  the  P  -  f>  method  was  made  by 
Dunnl3  for  plate  studies.  His  analytical  derivation  was  simplicity 
itself. 


He  wrote  a  functional  relationship  for  P,  the  applied  load,  and  for  the 
displacement  parameter  fo/  X.  He  argued  that  since  P  is  independent  of 
the  direction  in  which  the  sheet  buckles,  P  must  be  an  even  function  of 
fo/\.  Thus,  we  may  write  a  Taylor' series : 


P"  ,  .2  P"  1  ' 

P  =  P  +  (f  A)  + 

O  £7-  O/K  TJ, 


<f=A> 


(7) 


f  pi  I  pi  I  I  I 

Hence,  if  o/X  is  defined  as  u  and  the  coefficients  _  ,  ,  etc., 

are  defined  as  A,  B,  C,  etc.,  2'  P 


P  =  P  +  Au  +  Bu  + 
o 


(8) 


Thus,  for  small  values  of  u,  the  resulting  curve  is  very  nearly  a 
straight  line  and  the  intercept  with  the  load  axis  should  be  the  critical 
load.3 


Dunn  applied  his  method  to  the  analysis  of  data  which  he  had  obtained  from 
compression  tests  on  flat  panels  restrained  by  torsionally  weak 
stiffeners.  An  example  of  his  results  is  shown  in  Figure  11.  No 
correlation  was  made  with  the  Southwell  process.  It  was,  in  fact,  common 
contention  at  that  time  that  the  latter  method  is  inapplicable. 

13  Ik 

Following  the  research  of  Dunn,  Farrar  made  use  of  the  method  for 
interpreting  test  data  on  plates  under  axial  compression.  He  was 
concerned  with  plates  both  simply  supported  and  restrained.  Ifce  next 
application  appears  to  have  been  that  of  Hunk,  as  referenced  by  Hemp  and 
Griffin.  5  So  far  as  we  are  aware,  no  other  application  of  the  method 
was  made  until  recently. 

It  is  interesting  to  note,  however,  that  Donnell,  ^  in  his  classic  paper 
on  the  application  of  the  Southwell  method  performed  an  analysis  from 
which  both  the  o  and  the  b/P  versus  6  methods  can  be  obtained.  He 
showed  by  the  application  of  large  displacement  analysis  to  the  square 
plate  in  compression  that  the  relationship  between  the  displacement 
functions  is 

1  +  ^  (6  +  26p)(6  4  6Q)i  (9) 


cr  6+6 


8t 


\ 
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If  the  stipulation  is  made  that  _6  and  _o  should  be  small,  then  this 
equation  reduces  to  t  t 


Ik 


0 


7 - 1 - 1 

5  10  15  20 

P/A  X  10“3  (lb/sq  in.) 


Figure  11.  Dunn's  Experimental  Curves  for  Deter¬ 
mining  the  Buckling  Stress  of  the  Sheet. 
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(10) 


How  equation  (10)  is  clearly  identical  to  equation  (l),  the  small  die- 
placement  equation.  When  the  condition  that  6  /6  be  small  is  imposed, 
equation  (9)  simplifies  in  a  different  fashion0 to  become 


P  =  P  (1  +  y62) 
cr 


(11) 


where 


(12) 


This  is  the  large  displacement  formulation.  It  is  identical  in  character 
to  the  parabolic  approximation  to  the  elastica. 


The  validity  of  the  two  relationships  for  the  panel  in  axial  compression 
was  demonstrated  in  a  series  of  tests  made  on  square  fiber  glass  resin 
panels.  A  typical  load-displacemgnt  curve  is  given  in  Figure  12.  When 
the  appropriate  6/_  -  6  and  P  -  l>  curves  are  constructed.  Figures  13 
and  14,  the  linear^  relationship  is  clearly  seen  for  both  cases.  Moreover, 
the  values  of  P  determined  from  the  slope  and  intercept,  respectively, 
are  in  good  agr££ment. 


The  influence  of  initial  imperfection  is  readily  demonstrable  for  this 
case.  A  panel  similar  to  that  previously  described  was  tested  with  a 
normal  force  applied  at  its  center.  The  several  lgad-displacement  curves 
are  portrayed  in  Figure  15  •  The  6/P  -  6  and  P  -  6  curves  which  corres¬ 
pond  to  these  load  displacements  are  given  in  Figures  l6  and  17.  We  note 
that  the  Southwell  lines  are  parallel,  indicating  that  the  critical  load 
is  not  influenced  by  the  imperfection.  The  lines  on  the  P  -  6  plot  are 
also  seen  to  Intersect  on  the  F  axis,  indicating  that  the  critical  loads 
determined  by  the  large  displacement  methoa  are  the  same.  A  comparison 
of  the  values  determined  by  the  two  methods  indicates  that  they  are  in 
excellent  agreement.* 

It  has  already  been_shown  experimentally  that  the  Southwell  Plot  applies 
to  plates  in  shear.  The  applicability  of  the  large  displacement  process 
will  now  be  demonstrated  analytically  and  experimentally. 


*In  the  experiments,  bending  strains  were  used  instead  of  normal 
displacements.  In  this  class  of  problem,  since  the  units  of  the 
displacement  measurement  are  Irrelevant,  the  two  parameters  are  in  1:1 
correspondence. 
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17 


Figure  12.  Axial  Load  Versus  Bending  Strain 


20 


Figure  15.  Load  Versus  Bending  Strain  for  Tests  on  a  Fiber 

Glass  Panel  Under  Uniaxial  Compression  and  Lateral 
Load  Applied  at  Center. 


Figure  16.  Southwell  Plot  for  Data  of  Figure  15  • 
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The  analysis  is  baaed  upon  the  Kansan* Donna  11  finite  dsflactlon  theory 
for  which  the  compatibility  equation  ia 


r|s  Y  ■  VVJ 


K-  *  **  -  (  1  ♦  * 


for  the  imperfect  plate.  Ihe  assumption  will  be  made  that  the  elaetlc 
lateral  deflection  due  to  load  and  the  initial  iz^erfrction  shape  are  given 
hy 


W  =  yyA  Bin  512  .in 
LL  L  b 

W-TTa*  .in5E.inSM 

Li  Li  ISD  L  0 


Thus,  the  internal  strain  energy  U  can  be  expressed  ac 

_b  ju  r  _  1 


u '  I A  d  ( ^ ' 


for  the  singly  supported  plate. 

It  can  be  shown  from  equations  (13/  through  (13)  tha  equilibrium  in  the 
place  of  the  plate  is  satisfied  when 


i*. 

8L  b 


m  n  \  a 


2mrx  b 

cos  -  ♦  -»  cos 


2aiyr 

b 


Thus,  from  equations  (lU)  through  (17),  it  follora  that 


n  ♦  z]k  <*»/<**• **•«>* 


The  potential  energy  V  of  the  shear  If  is 

v  t 


v-i 


vff 


(w  ♦  W),X(W  ♦  w'),y4xdy 


By  substituting  equations  (l4)  and  (15)  into  equation  (19)  and  carrying 
out  the  Integrations,  ire  get 

T  4  A Wtt  ^5  r 

a  npq  (p  -  ■  )  (n  -  q  )  (2C 


where  m,n,  p,  and  q  are  Integers  such  that  map  and  n  t  q  are  odd  uuaberu. 
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When  the  normal  variational  principle  is  applied,  ve  deduce  after  appro¬ 
priate  algebraic  manipulation  that 
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p  q 

From  equation  (2l),  it  can  be  shown  that 

...  2n  ,4_4 

N  =  - 22—  (N  ) 

xy 


(21) 


4  4 


(V  +  2Ai 

(»V  ;  A2)2 


(22) 


The  experimental  work  to  verify  this  theory  is  taken  from  a  study  made 
by  Gerrard}-®  The  shear  stress  versus  bending  strain  diagram  is  shown  in 
Figure  18.  The  corresponding  small  and  large  displacement  plots  are 
given  in  Figures  19  and  20.  The  critical  stresses  derived  from  these 
curves  are  in  the  ratio  of  1:0.988,  being  4.66  ksi  and  4.55  ksi,  res¬ 
pectively.  This  experimental  value  is  93  percent  of  the  theoretical  value. 
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THERMAL  BUCKLING  OF  PLATES 


The  method  can  also  be  applied  in  the  case  of  thermal  buckling  of  a 
circular  plate  with  a  central  hot  spot.  This  question  was  first  discussed 
by  viueinec  19, whose  analytical  results  and  confirmatory  test  data  are  given 
here.  In  the  analytical  study,  ^ueinec  assumed  the  temperature  distri¬ 
bution  and  the  edge  conditions  were  perfectly  axis-symnetric .  In  the 
initial  condition,  the  plate  was  considered  at  uniform  temperature  and  free 
from  initial  stress.  The  effects  of  gravity  were  ignored  as  was  the 
variation  in  temperature  through  the  thickness  of  the  plate. 


He  concluded  from  his  analysis  that  if  the  displacements  were  small  and 
the  initial  imperfection  shape  and  the  final  deformation  mode  were  similar, 
the  displacement  at  the  center  could  be  related  to  the  initial  inperfection 
anplitude  by  the  equation 


W. 


w  = 


lo 


crit 


(23) 


This  expression  is  very  similar  to  the  Southwell  equation  for  the  strut. 
On  the  other  hand,  if  the  initial  and  buckled  shapes  were  different,  then 
the  relationship  is 


crit 


=  1  + 


1C  ^ 

h  *  w 


(24) 


where  is  a  coefficient  whose  values  depend  on  the  initial  shape. 

The  analysis  for  large  deflections  showed  that  in  this  case  the  deflection 
at  the  center  of  the  plate  and  the  tenperature  at  that  point  are  related 
to  the  critical  tenperature  and  a  constant,  Kg >  (which  is  associated  with 
the  deformation  mode)  by  the  expression 

To  -  Tcrlt  _/Vf 

Tcrit  V  )  (25) 

It  is  seen  from  this  formula  that  since  the  Tcrj_t  is  independent  of  Young's 
modulus  the  plate  deflections  are  also  independent  of  Young's  modulus  and 
depend  only  on  the  Poisson  ratio,  p,  of  the  material. 

Clearly,  when  the  experimental  data  giving  W0  as  c  function  of  T0  are 
available,  these  equations  can  be  used  to  determine  Tcrit» 
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LATERAL  INSTABILITY  OF  DEEP  BEAMS 


20 

In  recent  research,  Way  has  demonstrated  the  applicability  of  both 
large  and  small  displacement  methods  to  the  problem  of  lateral  instabi¬ 
lity  of  a  deep  beam  loaded  by  concentrated  force  lying  in  the  plane  of 
the  web  and  passing  through  the  centroid.  He  obtained  excellent  agree¬ 
ment  between  the  two  methods.  The  load  displacement  which  he  determined 
is  given  in  Figure  21.  The  6/p  versus  6  plot  is  shown  in  Figure  22,  and 
the  P-62  is  shown  in  Figure  23.  The  critical  load,  levels  determined 
from  these  curves  are  26l  grams  and  253  grams  which  compare  well  with 
the  theoretical  value  of  260  grams . 
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Figure  23.  P  Versus  6  for  Data  of  Figure  21 


CONCLUSION 


This  technique  gives  identical  results  to  those  obtained  from  the  well- 
established  Southwell  process  but  is  easier  to  apply  in  many  cases.  Both 
methods  have  their  own  peculiar  restrictions.  In  the  large  deflection 
case,  the  load  levels  must  be  in  excess  of  the  theoretical  critical  value; 
but  it  is  clear  from  the  data  available  that  small  amounts  of  local  yield¬ 
ing  do  not  seriously  influence  the  problem.  For  the  Southwell  method  to 
apply,  the  motions  must  be  small  and  the  initial  imperfection  likewise 
restricted  for  many  cases.  Thus,  for  the  small  displacement  approach, 
there  is  an  implied  requirement  for  a  relatively  good  quality  test  vehicle, 
a  requirement  of  considerable  less  importance  for  large  displacement 
considerations.  However,  in  all  cases,  the  user  must  be  cautioned  against 
the  application  of  either  method  when  the  data  are  such  that  curve  fitting 
techniques  would  be  needed. 
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of  experimental  results  on  a  wide  range  of  simple  structures.  The  first  of 
these  applies  when  the  deformations  are  small!  NIt  is  that  the  Initial  «.nH 
subsequent  deformations  are  related  to  the  actual  and  classical  instability 
loads  by  the  expression 


It  is  likewise  shown  that  when  the  deflections  are  large  the  loads  are 
in  excess  of  the  instability  load>!’the  relationship  is  modified  to 

P  =  Pj(l  +  Y62) 
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